Heterogeneities in systems with quenched disorder 
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We study the strong role played by structural (quenched) heterogeneities on static and dynamic properties of 
the Frustrated Ising Lattice Gas in two dimensions, already in the liquid phase. Differently from the dynamical 
heterogeneities observed in other glass models in this case they may have infinite lifetime and be spatially 
pinned by the quenched disorder. We consider a measure of local frustrationshow how it induces the appearance 
of spatial heterogeneities and how this reflects in the observed behavior of equilibrium density distributions and 
dynamic correlation functions. 



I. INTRODUCTION 

Much of the work in the theory of the glass transition has 
been concentrated in the precursor phenomena present in the 
high temperature (or low density) phase on approaching the 
glass transition. The (equilibrium) dynamics in this region 
shows two qualitatively different regimes: a short time re- 
laxation associated with the rattling of particles inside cages 
formed by their neighbors and a long time structural relaxation 
which, as density increases (or T decreases), takes longer and 
longer times |1,2]. In supercooled liquids, the presence of 
a crossover temperature below which exponentially decay- 
ing correlation functions typical of the liquid phase become 
stretched has been observed f^. The origin of this stretch- 
ing of relaxations is associated to the gradual appearance of 
heterogeneities in the space/time domain. Understanding the 
emergence of these heterogeneities is at the heart of the na- 
ture of the glass transition and has beenjhe source of a lot 
of theoretical and experimental work |4, 5, 6]. Neverthe- 
less, a unifying microscopic theory is still lacking and only a 
few universal properties related to heterogeneities have been 
found. 

In real space, dynamical heterogeneities coiTespond to par- 
ticles with different mobilities. In computer simulations of 
simple glass formers the presence of fast and slow particles 
has been observed. The typical time scales of the two types of 
particles can be identified for example by examining the be- 
havior of the non-gaussian parameter associated with the van 
Hove correlation function |7, 8]. This leads to the identifi- 
cation of a growing time scale which seems to diverge at the 
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mode coupling transition temperature (Tmct)- This charac- 
teristic time arises naturally when computing a suitable time 
dependent four point correlation function or dynamical non- 
linear susceptibility As the temperature approaches 
Tmct from above, this dynamical response shows a growing 
peak at a characteristic time t* and then relaxes at long times 
to its equilibrium value. This growing characteristic time of 
the four point functions has been associated with the presence 
of long lived dynamical heterogeneities (clusters of very slow 
or immobile particles). 

From a landscape point of view, the emergence of stretch- 
ing seems to be consequence of the gradual confinement of the 
system in low (free) energy regions of phase space surrounded 
by growing barriers as temperature decreases. The confine- 
ment is due in part to the dramatic decrease of escape direc- 
tions in configuration space on approaching the glass tran- 
sition 1 1 L .12. .13il . In simple models of binary mixtures it 
has been found that stretched relaxations are associated with 
long lived superstructures or "metabasins" in configuration 
space 1 14]. The relaxation from single metabasins shows a 
stretched behavior and different metabasins relax with differ- 
ent characteristic times, giving rise to the overall stretching 
observed for example in the a region of the incoherent scat- 
tering function. 

At variance with the enormous amount of work on the prob- 
lem of heterogeneities in glasses, much less is known about 
them in systems with quenched disorder fisl ITS [TtII . Dif- 
ferently from true glasses, where the inhomogeneities are dy- 
namically created by the self induced local frustration and per- 
sist during a given timescale, the quenched disorder induces 
structural heterogeneities that may have infinite lifetime. In 
this sense they can be called quenched heterogeneities. Since 
the inhomogeneities are pinned by the disorder, they may be 
easier to detect and characterize than in systems with self in- 
duced disorder. Thus, the posed question is to what extent 
these heterogeneities influence the relaxation properties and 
how different regions of the sample differ from each other. 
Moreover, how do they affect the different degrees of freedom 
present in the system (like particles and spins)? 
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To answer these questions, in this work we address the 
problem of heterogeneities in systems with quenched disorder 
by focusing the analysis on the behavior of a two dimensional 
version of the Frustrated Ising Lattice Gas (FILG) which has 
been proposed as a simple lattice glass model [18]. Neverthe- 
less it is known that the presence of quenched disorder in the 
definition of the model introduces some features not typical 
of real glasses. In particular it acts as a pinning field for het- 
erogeneities, some of which can have infinite lifetime. This 
is reflected in the behavior of dynamical correlation functions 
and in the non-linear susceptibility and compressibility as will 
be shown below. In this context we will consider a measure 
of local frustration lITslflTil . show how it induces the appear- 
ance of spatial heterogeneities and how this reflects in the ob- 
served behavior of equilibrium density distributions and dy- 
namic correlation functions. For the sake of comparison with 
other glassy systems without quenched disorder we also dis- 
cuss the behavior of the non-linear susceptibility associated 
with the spin degrees of freedom and of the non-linear com- 
pressibility of the density variables. Finally, an analysis of the 
universal behavior in terms of hole variables recently intro- 
duced 1 19, 20] shows that this system seems to be in a dif- 
ferent universality class than the previously investigated glass 
models without disorder. 



n. THE MODEL 

The Frustrated Ising Lattice Gas jll|2ll|23|2l is defined 
by the Hamiltonian 

i7 = - J ^ {eijaiffj - l)ninj - ^^n^ . (1) 

<ij> i 

At each site of the lattice there are two different dynami- 
cal variables: local density (occupation) variables ni = 0,1 
(i = 1 . . . N) and internal degrees of freedom, ai — ±1. The 
usually complex spatial structure of the molecules of glass 
forming liquids, which can assume several spatial orienta- 
tions, is in part responsible for the geometric constraints on 
their mobility. Here we consider the simplest case of two pos- 
sible orientations, and the steric effects imposed on a particle 
by its neighbors are felt as restrictions on its orientation due to 
the quenched variables connecting nearest neighbors sites. 
The first term of the Hamiltonian ensures that when J is large, 
any two neighboring particles will have their spins satisfying 
the connection between them. Finally /i represents a chemical 
potential ruling the system density (at fixed volume). 

This model has been studied mainly through simula- 
tions J 18,^, '23'] in 3d and analytically in mean field ]24,'23]. 
In the low density (high temperature) phase, the system behav- 
ior is liquid like, the dynamics is fast, time-translationally in- 
variant and obeys the fluctuation-dissipation theorem. By in- 
creasing the density (or lowering the temperature), the system 
has a spin glass transition characterized by a divergence of the 
static nonlinear susceptibility. This transition point is located 
quite close to the point where the system suffers a dynami- 
cal arrest, where the particles are no longer able to diffuse. 
After a sub-critical quench, the dynamics is slow and history 



dependent. Besides that, it presents aging and violates the 
fluctuation-dissipation theorem (we refer the interested reader 
to ]23] and references therein). 

III. STATIC PROPERTIES 

Since bonds are created at random (with probability 1/2 of 
being ±1), the system is not homogeneous at the lenghtscale 
of the lattice spacing. Each site belongs to several minimal 
plaquets and there is a probability that all (or neither) of these 
are frustrated. This would influence the dynamics at a lo- 
cal level since the system is no longer invariant under spatial 
translations. Thus we may ask which are the consequences of 
these quenched heterogeneities on the equilibrium properties 
of the system. 

When evaluating the partition function at finite density, the 
configurations in which particles close any frustrated loop will 
have zero probability. The average density will no longer be 
site independent and the density distribution P{p) will be in- 
homogeneous, i.e. will show spatial fluctuations. Homogene- 
ity will only be recovered when p ^ 0. For small densities, 
the distribution will be gaussian while deviations from gaus- 
sianity should be evident as the density increases. The forbid- 
den configurations introduce spatial inhomogeneities that are 
long lived as they derive from the quenched random underly- 
ing connections. In other lattice models, as well as in struc- 
tural glasses, the disorder is not quenched, what restores trans- 
lational invariance if the measure time is much larger than the 
dynamical heterogeneities characteristic time. An example of 
a system with quenched but not random frustration, that does 
not present these inhomogeneities is the fully frustrated ver- 
sion of the FILG. 

To quantify this we use the (local) coarse-grained quantity 

^■'^E n (2) 

where V{i) are the minimal plaquets containing the site i. 
That is, for every site, we consider the 4 minimal plaquets 
that contain that site (12 in d = 3) and we check how many 
of these plaquets are frustrated, i.e. the product of the connec- 
tions in the loop is negative. This local frustration parameter 
was originally considered in the context of dynamical hetero- 
geneities in refs. fTsl flTil . In this nearest neighbor plaquet 
approximation sites that are maximally frustrated have all of 
its 4 plaquets with negative products and sites with no frustra- 
tion have all plaquets with positive products. In other words, 
7i may assume values between —4 and 4. Figure^Ji shows a 
configuration with inhomogeneous distribution of local frus- 
tration. 

We now consider how these structural clumps change the 
system's behavior In fig. [Q) we show, for the same realiza- 
tion of disorder of fig.^, only the sites that have an averaged 
occupation greater than 0.9. In this figure, fi — A, and the 
total average density is roughly 0.74, well below 0.9. The 
presence of preferred positions is a direct consequence of the 
existence of islands of low frustration, as can be seen by com- 
paring figures^ and^. As the chemical potential increases. 
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FIG. 1: a) Example of a {7^} configuration. Black sites are those 
whose four neighbor plaquets are non frustrated (7^ = 4) while white 
sites have all frustrated (7^ = — 4). Intermediate values of 7^ are rep- 
resented by different degrees of grey, b) Sites whose average density 
is higher than 0.9 for n — 4, corresponding to a global p ~ 0.74. 

the number of these sites also increase, as expected. These 
heterogeneities make the system not translationally invariant 
even in the low density liquid phase. This can be made clear 
by measuring the distribution P{pi) of local averaged densi- 
ties Pi — (rii) as shown in fig.|2] In homogeneous systems 
all sites have the same average density and the distribution, 
for a finite number of measurements, is a gaussian centered in 
this value. Examples might be a non disordered version of the 
model considered here (with fully frustrated connections) or 
models for structural glasses without quenched variables. In 
the last case, the annealed character of the frustration allows 
the homogeneity to be recovered again. The broadness of the 
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FIG. 2: Distribution of local densities for several values of p. and 
L = 20 in 2d. Notice the multi-peak structure that becomes more 
apparent as p increases. 

curves does not seem related neither to finite average times 
nor to finite lattice sizes, but rather to the fact that, when not 
limited to the nearest neighbor plaquet approximation, there 
are several values that the local frustration may assume gen- 
erating in this way the broad distributions seen in the above 
figures. 

To make clear the origin of this broad distribution, we 
considered the density distribution P{pi,ji) separately for 
each value of the local frustration 7^, which has a well de- 
fined density around which its particular distribution is cen- 
tered. The result is seen in fig. |3l Moreover, the separation 



between these curves increases as the density increases and 
the deviation from gaussianity becomes apparent. We may 
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FIG. 3: Distribution of local densities for p — 3 and L = 20 in 
2d. The bold curve is the sum of the several others, P{pi)- From left 
to right, we have 4, 3,. . . ,0 locally frustrated loops. The distribution 
has an internal structure. Inset: The distance A between the average 
densities for extreme values of 7, P{pi, 4) and P{pi, —4). 

define A as the difference between the average densities of 
the peaks associated with the lowest and highest values of 
7, A = (pi)^ - {pi)_^ where the averages (...)^_ are over 
P{Piili)- When there is no difference, A = 0, the parti- 
cles do not feel the underlying landscape and all sites have the 
same average equilibrium density. In this model this happens 
for very low densities. As the density increases, the dynamics 
becomes landscape influenced and P{pi) has a more complex 
structure, no longer being a simple gaussian. This is exempli- 
fied in fig.l^for p = 3. One can also note that the distribu- 
tions are asymmetric around % = 0: for example, in the case 
of P{pi, —2) and P{pi, 2) the former is broader and shorter 
while the latter is higher and more concentrated. A possible 
explanation is that sites with negative values of ji are more in- 
fluenced by higher order plaquets than sites with positive 7^, 
thus having a broader distribution. 

In d = 3 the overall picture is similar to that of d = 2. 
Nevertheless some differences are evident when comparing 
figures |2] and 13 and one may wonder if the minimum plaquet 
approximation used in the definition of the local frustration 
parameter works in d = 3 as it does in d = 2. Still within 
the same approximation, the distribution of local densities for 
different local frustrations presents also different peaks (not 
shown), a manifestation of the lack of translational invariance 
also in d = 3. 

Recently, a description of the jamming transition has been 
introduced that unveiled some universal mechanism leading 
to the dynamical arrest that happens at pc- The notion of 
hole is introduced: empty sites that have at least one neigh- 
bor particle able, due to the energetic or to the kinetic con- 
straints, to jump to the initial empty site. We measure, for 
several^xet/ densities |33], the density of holes, i^. It has been 
conjectured, and supported by numerical simulations on some 
lattice models 1 19, 20], that the diffusivity depends on i' as 
D ^ {v — vq)^ where vq is a possible residual density of 
holes at pc (rattlers) that do not contribute to the diffusion. 
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FIG. 4: Distribution of local densities for several values of /i and 
L = 10 in 3d. Differently from the 2d case, the curves do not seem 
to present internal structure. 



We show our results for the FILG in d — 2 in fig. 13 We 
obtain D ^ {v — O.OSG)^ "*. The Kob-Andersen is another 
example of model where vq ^ 0. Some things are remark- 
able: first, we notice that the FILG belongs to a different uni- 
versality class. Second, the value of v at pc is much larger 
than in other models. Notice also that the behavior of v near 
Pc is linear. The fact that the FILG does not belong to the 
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FIG. 5: The fraction of holes i/ as a. function of p — pc for L — 50. 
Notice that u does not vanish at pc ~ 0.79. The behavior near pc is 
linear with the density. Inset: diffusion coefficient as a function of 
pc — p. The diffusivity goes to zero as (pc — p)'^ with ~ 2.4. 

same universality class of several other models is very inter- 
esting and the reason may be related to the inhomogeneity of 
the underlying connections structure. Indeed, we measure the 
density of holes accordingly to the ji of the hole site. For 
small values of p there is no difference between these values 
and all curves collapse. The fact that for higher densities these 
curves fall apart is an indication that the holes also have pre- 
ferred sites and tend to form small clusters. Thus, near pc, 
the mechanism leading to diffusion is no longer the pairwise 
collisions between holes but a more complicated process in- 
volving these small clusters, what may lead to a failure of the 
i/-square law l34ll . Further work along this line has to be done 
in order to clarify this issue. 



IV. DYNAMICAL PROPERTIES 

In Ref.^3] the nonequilibrium behavior of the model was 
shown to present interesting aging properties similar to those 
observed in glass forming systems. Here we address equi- 
librium density correlations instead and the role of hetero- 
geneities in the low density phase. The density autocorrela- 
tions are defined by 



Cnit) = 



n^{t)n,{0)) - {n,{0)y 



(3) 



and Cn{t) = c„(t)/c„(0), where the averages are both over 
samples and initial states. The densities are also averaged 
over the thermal history in each sample. Notice that only in 
the case where the system is homogeneous, that is, the av- 
erage quantities are invariant under spatial translations, we 
can substitute the last term in the numerator by p^, other- 
wise i^'i)'^ 7^ "j)^- In the region of densities consid- 
ered here, these coiTelations do not present the usual two step 
dynamics observed in other models as can be seen in fig. |6l 
Nevertheless two regimes can be clearly observed and all the 
curves can be fitted by stretched exponentials (not shown) in 
the first regime. The second regime at longer times observed 
for the largest densities can also be fitted by stretched expo- 
nentials with other parameters and probably will tend to de- 
velop a plateau at higher densities. 




FIG. 6: Equilibrium particle correlations C„{t) for several values 
of p for a dynamics with spin flip, particle creation-destruction and 
diffusion. 

For the same values of the chemical potential or mean den- 
sity shown in fig.|6j the diluted spin correlations show a dif- 
ferent behavior The correlations between spins on occupied 
sites are defined as 

Cnsit) - <"^W'^^W"^W5.(0)) . (4) 

i 

The results for several values of p are shown in fig.Qtogether 
with the best fits to stretched exponential decays for the long 
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time regime. Interestingly, the correlations C„ and Cns de- 
cay with very different timescales: density correlations de- 
cay much faster than the diluted spin correlations. This can 
be understood as follows. Particles can be created/destroyed 
and diffuse through the lattice, forming clusters of neighbor- 
ing particles. Once they belong to the same cluster, the nature 
of the interaction enforces them to have spins satisfying all the 
bonds in the cluster, otherwise they should move apart. Thus, 
when a site contributes for the correlation at different times, 
the contribution will be positive, unless there is enough time 
for the whole cluster to flip all its spins. Moreover, the perco- 
lation transition is just below /i = 3 and the considered clus- 
ters are large, making the decorrelation time quite large. The 
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FIG. 7: Equilibrium diluted spin correlations C„s{t) for several val- 
ues of fi for a dynamics with spin flip, particle creation-destruction 
and diffusion. The symbols are the same as in fig.|6| 

diluted spin correlations start to develop an incipient plateau 
as the value of the chemical potential increases. We expect 
that by further increasing the value of /i, this plateau will be 
more noticeable and the characteristic decay time will tend 
to diverge at the dynamical transition. In fig.|S]we show an 
Arrhenius plot of the relaxation times showing that they tend 
to diverge for /i ^ oo. This is different from the behavior 
observed in the three dimensional model |26], which shows a 
power law divergence at a finite value of fj,. Thus, while the 
three dimensional model presents a spin glass transition at a 
finite chemical potential the two dimensional model only has 
a transition at /i ^ oo. 

Both correlations can also be measured only for sites pre- 
senting a given value of the frustration parameter 7;. The 
curves for the 5 possible values (in d = 2) are shown in figs.|9l 
and^Jfor /i = 3. Fits to stretched exponential decay are also 
shown for the diluted spin correlations. They are normalized 
to unity but, recalling fig.|3 the number of sites contributing 
to each one of the curves is different. The behavior for C„ (i) 
is very different and almost all curves give the same corre- 
lation, only sites with 7^ = 4 differ slightly from the others 
(and their contribution to the overall curve is small). On the 
other hand, for Cns (t) the curves are strongly resolved and we 
notice that the dynamics is faster in those sites in highly frus- 




FIG. 8: Relaxation times of the diluted spin correlations as a function 
of the chemical potential in a semi-log plot. 



trated regions. This may be due to the fact that since particle 
creation-destruction is allowed, particles appearing in these 
regions have a higher probability to go elsewhere or be de- 
stroyed (as these are regions of less than average densities), 
what makes the decorrelation time small. On the other hand, 
low frustration regions are those preferred by the particles and 
correlations take much more time to decay (as these are re- 
gions of larger than average densities). The total correlation 
(not shown) follows closely the 7^ = or mean frustration 
curve. Note that all local correlations show stretched exponen- 
tial relaxations irrespective of the value of the local frustration 
being high or low. In this scenario the stretched relaxation 
does not emerge as a consequence of some kind of convolu- 
tion of local exponential relaxations with different time scales. 
At this level of description, the stretching, although heteroge- 
neous due to different local frustration, is intrinsic. As the 
density increases this scenario does not change qualitatively. 
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FIG. 9: Equilibrium particle correlations Cn{t) for the five possible 
values of ji and /i = 3. The line is the full correlation, from fig.|6| 
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FIG. 10: Equilibrium dilute spin correlations Cns (t) for the five pos- 
sible values of 7^ and jj. — 3. The symbols are the same as in fig.|9| 
Notice that low frustration sites already show a small plateau, al- 
though in the full correlation (fig.0, this is not yet noticeable. 



FIG. 11: x'ii't) as a function of time for different chemical poten- 
tials. For increasing chemical potential (hence densities) the peak 
increases and shifts to longer times. 



It is interesting to remark that for high values of /i the decay 
of the diluted spin correlation starts to develop a plateau and 
the sites responsible for this behavior are those that are less 
frustrated. The effect is not much pronounced because the 
summed up effect of 7^ = —2 and sites is dominant. 

As the dynamical heterogeneities have attracted much at- 
tention recently, several quantities have been devised to quan- 
tify them Pd . Among them, of particular interest is the four 
point correlation (or dynamical nonlinear response) defined as 

xm ^ N [{Cl{t)) - {QMf) • (5) 

Away from the transition, this function presents a maximum at 
a time t* that is a measure of the timescale during which the 
particles forming the heterogeneity are coiTelated. The long 
time behavior of this quantity is a measure of its equilibrium 
value. In figure[^we plot X4 (0- Analogously to what hap- 
pens in other glass models it presents a broad peak that shifts 
to longer times and gets higher as the density increases. Nev- 
ertheless an important difference from other models is that af- 
ter the peak the non-linear compressibility does not decay too 
much and develops a density-dependent plateau at long times. 
The persistence of density correlations is a clear signature of 
the presence of very long lived heterogeneities, a consequence 
of the pinning effect of the quenched disorder present in the 
model. The possible divergence of the equilibrium (infinite 
time) limit of the non-linear compressibility should point to 
the existence of a thermodynamic transition associated with 
the density degrees of freedom not yet observed in this model. 

Analogously, for diluted spin correlations we can evaluate 
eq.|5]using Cns (t) in the place of C„ (t). The resultant X"* {t) 
is seen in fig. [21 The behavior is similar to that observed 
for the density variables, but instead of a peak we observe a 
crossover time after which this non-linear susceptibility satu- 
rates in a finite density dependent value. In this case the height 
of the plateau seems to diverge much more rapidly than for the 



compressibility giving more plausibility to the presence of a 
thermodynamic transition in the spin variables, probably for 
/i ^ 00. The crossover to a plateau, instead of the presence 
of a peak, clearly points to the appearance of a persistent cor- 
relation length associated with the spin degrees of freedom. It 
would be interesting to check if these persistence is a manifes- 
tation of a growing structural correlation length which might 
diverge at an equilibrium spin glass transition similar to that 
present in d=3. 
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FIG. 12: x^'it) as a function of time for different chemical poten- 
tials. Note here the saturation after a crossover time and the stronger 
growth of the subsequent plateau at variance with the behavior of 

xm in fig. mi 



V. CONCLUSIONS 

In this paper we considered the role played by hetero- 
geneities present in a model with quenched disorder. Dif- 



7 



ferently from self-induced frustration models here the hetero- 
geneities may be pinned by the local frustration and this al- 
lows a simple characterization of them. The effects of these 
pinned heterogeneities on time correlation functions and dis- 
tribution of local densities have been analyzed. The system 
forms clumps in the Uquid phase even when the particles ei- 
ther do not interact or strongly repel each other f27"]. More- 
over, these clumps are long-lived and it would be interesting to 
study their size distribution and how they are correlated f28ll . 

A simple measure of local frustration f]3\ was considered 
and a strong correlation with the inhomogeneous density dis- 
tribution was found. The use of 7^ allowed to decompose the 
local density distributions in components related to the degree 
of local frustration. However, the use of 7; as a measure of the 
local frustration considering only the minimal plaquets has to 
be considered as a first order approximation. Indeed, higher 
order plaquets, involving larger loops, have to be taken into 
account in order to obtain the full dependence of pi on 7^. 
For example, in the d = 2 case, there exists 12 plaquets of 
second order (those having 6 bonds). It is an open problem 
to determine to what extent the present results may change 
by considering the effects of larger loops in 7^ and how they 
determine the average local density and affect the relaxation 
properties. It was shown also that dynamic correlations can 
be decomposed accordingly to the local values of frustration 
and in each case the relaxations are essentially stretched al- 
though the degree of stretching is of course dependent on the 
degree of local frustration. It would be interesting to study 
the possible relation between local correlation functions in the 
out-of-equilibrium regime and the violation of the fluctuation- 
dissipation theorem Il6l 

m, 113, EH 

with the local frustration 

as defined by the parameter 7^ . 

Further insight into the relation of the FILG with other 
models of glasses and spin glasses can be gained by analyz- 
ing non-linear responses. We showed that two such responses 



associated with different degrees of freedom in the model be- 
have differently in some important aspects. At variance with 
models of glasses, the spin non-linear susceptibility shows a 
crossover time after which the response saturates instead of 
the characteristic peak observed for example in models with 
constrained dynamics. Also the rapid growth of the plateau at 
long times suggests the developing of a structural correlation 
length associated with a possible spin glass transition. This 
remains to be confirmed. The non-linear compressibility in- 
stead shows a broad peak but with a modest further relaxation 
to a density dependent plateau. In this case one can associate 
this plateau with the formation of very long lived quenched 
heterogeneities. From the observed behavior of these non- 
linear responses one can conclude that the FILG does not be- 
have neither as a proper spin glass nor as a glass sharing some 
characteristics with both kinds of systems. 

Although the model considered here has been originally in- 
troduced as an attempt to obtain a finite dimensional lattice 
system with glassy properties similar to structural glasses the 
main drawback is the lack of invariance under spatial trans- 
lations already in the liquid phase. There are several ways to 
restore this invariance. For example, by considering spatially 
coarse-grained variables or by including a temporal evolution 
of the bonds |21]. As it stands this model still may be in- 
teresting to study transport in disordered media (see |32] and 
references therein) which has been of much interest recently. 
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